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Abstract 

For the two dimensional randomly disordered three-state Potts ferromagnet, 

we find numerically that the critical exponent rj unchanges with the degree of 

disorder while those of the correlation length {y) and magnetic susceptibility 

(7) increase with it continuously. We discuss some consequences of the finding. 
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The effect of random quenched disorder to the critical behavior of a ferromagnetic system 
has been a subject of intensive studies. According to a heuristic argument by Harris 0, that 
has been widely accepted over two decades, the disorder induces a new universality class 
only when the critical exponent of the specific heat (a p ) in the corresponding undisordered 
system is positive. 

A more rigorous approach f|| claims that in a D-dimensional disordered system obeys 

u>2/D. (1) 

The proof of Eq. ([[]), however, presumes the existence of so-called finite-size scaling correla- 
tion length whose relation to the intrinsic correlation length £ is by no means obvious, defined 
in terms of probability. Note that Eq.(|l|) combined with hyperscaling relation, a v = 2 — Du p 
(with Up denoting the value of u in the corresponding undisordered system), implies that 
for a disordered system with a p > u jumps discontinuously from the undisordered value 
Up. To be more specific, random disorder is supposed to induce a different scaling region 
near criticality while the rest of the scaling region is equivalent to that of the undisordered 
system (crossover from the undisordered to disordered). 

The 2D random bond disordered Ising ferromagnet is unique in that its continuum limit 
is found to be a certain tractable field theoretic model, i.e., O(N) Gross- Neveu model in the 
limit iV — ► 0. Accordingly it is predicted that the value of rj in this system does not vary 
with the strength of weak disorder 0, which has now been confirmed by various numerical 
methods even for strongly disordered case 0-f§. Albeit a p = 0, recent extensive numerical 
studies find that both u and 7 change with disorder, in disaccord with Harris criterion. 

Most numerical studies on the 3D random site diluted Ising ferromagnet, where a p ~ 0.1, 
report evidences for new universality class [^J,|TT],|T2j . Still is it not clear whether Eq.(|l|) is 



indeed satisfied or not ||11|| ; nor is it clarified whether critical exponents change continuously 
with the strength of disorder or jump to certain fixed values regardless of it. An accurate 
estimate of rj is far more difficult, for it is very sensitive to the choice of the critical point. 
Generally speaking, current understanding on generic randomly disordered ferromagnet is 



rather limited theoretically, numerically, and experimentally. 

In this Letter we report new finding of the extensive Monte Carlo studies on the 2D ran- 
dom bond disordered three-state Potts ferromagnet. The technical difficulty is significantly 
lightened in this system, since the exactly known value of a v is large enough (a p = 1/3) 
and the exact critical point can be calculated for certain realizations of disorder. Thus we 
obtain clean numerical evidence that the value of r\ unchanges with disorder while both v 
and 7 increase continuously. With the very same results found in the 2D disordered Ising 
ferromagnet ||, we conjecture that those features are generic in any randomly disordered 
ferromagnet s. 

The Hamiltonian of a g-state Potts ferromagnet with quenched random interaction 
(bond) can be written as 

H = - J] JxySa x a y , (2) 

<x,y> 

where the spin at site take on the values l,...,q (q = 3 in this work), 5 is the 

Kronecker delta function, and the sum is over all the nearest-neighbor pairs on the lattice 
system. The coupling J xy for each link < x, y > is selected randomly from two positive 
(ferromagnetic) values J and J' with probability p and 1 — p respectively. 
For p = 1/2 the system is self-dual with its self-dual point given by ]13| 



[exp(J)-l][exp( J') -!] = <? (3) 

We fix J = 1 and p = 1/2 without any loss of generality, and consider four different values 
of J', i.e., J' =1.0 (undisordered case), 0.9, 0.5, and 0.25. Accordingly, the inverse self- 
dual temperature (critical temperature) is given by p e = Ib(V3 + 1), 1.0588..., 21n(2), and 
1.8235..., respectively for J '=1.0, 0.9, 0.5, and 0.25. The values of the critical exponents 



for the undisordered case are exactly known |TJj], e.g., v v = 5/6 and 7 P = 13/9 (hence 
7^ = 4/15). 

Our raw-data for each J' are obtained by choosing a realization of random J', then 
running Monte Carlo simulations in Wolff's single cluster algorithm with periodic boundary 



conditions imposed on square lattice. Details of our single cluster algorithm for Potts model 
will be presented in a longer paper. For each realization, measurements were taken over 10 
000 configurations each of which is separated by 6-24 one cluster updatings according to 
autocorrelation time. The mean values over all the different realizations are our final data. 
The numbers of different realizations are typically 30-60, 150-200, and 250-350 respectively 
for J 1 = 0.9, 0.50 and 0.25. Larger number of realizations are required for smaller value of 
J' due to wilder fluctuations among different realizations. 

Our magnetic susceptibility and correlation length on a square lattice of linear size L are 
defined as 

XL = G(0), (4) 

& = ^m^ J ^ I) - 1 (5) 

where the two-point correlation function G(x) and its Fourier transformation G(k) are re- 
spectively defined as G(x) =< 5 aoax — l/q> and G(k) = Y^ x G{x) exp(kxi) with x\ denoting 
the first component of x. 

In order to determine the value of r\ we measured xl at each critical point (self-dual 
point) by varying L from 20 to 120. According to the standard theory of finite size scaling, 
Xl at criticality satisfies 

XL ~ L 2 ~\ (6) 

where r\ = 2 — 7 jv. Our results are summarized in Figure(l) showing undoubtedly that the 
value of 7/1/ remains unchanged irrespective of the value of J'. The x 2 fits yield ^=0.266(3), 
0.264(3), 0.267(8), 0.268(10), respectively for J' =1.0, 0.9, 0.5, and 0.25, to be compared 
with the exact value 0.266... (x 2 /N DF < 0.5 for all the J'). 

The value of 7 is determined by the x 2 fit of bulk x data measured at various values of 
the inverse temperature (3. The data with £((3) > 5 are considered in the fit to eliminate 
possible correction to the scaling 

X - t~\ (7) 
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where t = (/3 C — /3)//3 c . We used £l to monitor the finite-size effect in our measurements of 
bulk X] that is, the condition > 8 was always imposed for the extractions of our bulk 
data (see Figure(2)). The largest values of £ are roughly 20 ~ 30 for all the values of J', 
using L up to L=240. 

Our procedure for the undisordered case yields 7=1.442(3) in excellent agreement with 
the exact value 7=1.444... . The data for different J' are summerized in Figure(3). The 
value of 7 clearly increases with the strength of disorder: 7=1.44(1), 1.56(4), and 1.84(1) 
respectively for J'=0.9, 0.5, and 0.25 {x 2 /Ndf < 1)- Also we do not find any indication 
of crossover, as in the 2D randomly disordered Ising ferromagnets |||8| . In other word, the 
presence of random quenched disorder affects the entire scaling regime (Figure(3)). 

For the weakly disordered case (J 7 = 0.9) the value of 7 thus extracted is virtually 
the same value as in the undisordered case. We do not take this as an evidence that 
their values are strictly the same in those two systems; it rather appears that 7 increases 
indistinguishably mildly for a very weakly disordered system. 

Eq. (JT|) combined with our finding that 7/z/ remains a constant implies that 7 > 26/15 
in the two dimensional randomly disordered three-state Potts ferromagnets, in obvious dis- 
accord with our results for J' = 0.9 and 0.5. 

Our results contrast with some previous numerical results in the 3D disordered Ising 



ferromagnet JT0| fL5 |, but are consistent with the results in Ref. [12] where the critical 
exponent of the magnetization was found to vary continuously. The variant conclusions of 
the previous numerical studies on the 3D disordered Ising ferromagnet are perhaps due to 
some technical reasons: Without measuring correlation length, a priori, the finite size effect 
in the measurement of bulk data can hardly be monitored, nor can there be a criterion to tell 
whether or not a data point (/9,x(/3)) is in the scaling regime. Moreover, without knowing 
precise critical point one usually needs very broad range of the data for a reliable extraction 
of the parameters from a fit to Eq.(|7]). Our current investigation is free of all the drawbacks. 

Our results bear striking resemblance to the results of the 2D random-bond disordered 
Ising ferromagnet |§: There the data up to £ ~ 540 are analyzed for the similar values of 



J' to those here, and it is found that v and 7 for J' = 0.9 are virtually the same values as 
in the corresponding undisordered system while they are clearly increased for J' = 0.25 and 
0.1. It thus appears that the critical behavior of a disordered system is not characterized by 
the value of a p . 

Some time ago Suzuki proposed |]TB| the concept of weak universality to comprise the 
rigorously known critical exponents of the eight- vertex model into the concept of universality. 
According to it the value of v may change with the details of microscopic interaction whereas 
the value of 77 does not, so that a (weak) universality class is characterized by the value of 
1], but not by the value of v or 7. The presence of the random quenched disorder does 
not change the symmetry of the system although microscopic interaction changes. Thus, in 
view of the weak universality the value of 77 remains the same value as in the corresponding 
undisordered system whereas both v and 7 vary. 

Some non-trivial examples of the realization of the same weak universality class would 
be the D-dimensional XY and Ising models, for the Z(2) group in the Ising ferromagnet 
is the center group of the 0(2) symmetry in the XY model. It has recently been proved 
numerically that the value of 77 in the 2D XY model is indeed 1/4 jTT]] (as in the 2D Ising 
model), while r\ = 0.031(4) and 0.026(3) [0 respectively for the 3D XY and Ising models 
being equal to each other within the statistical errors. We conjecture that the similar 
pattern of weak universality should manifest itself even in a disordered system with a p < 0: 
the 3D disordered XY and Ising models should be in the same weak universality class as the 
corresponding undisordered models, albeit a p is negative (positive) in the 3D undisordered 
XY (Ising) models. 

To conclude: We have shown that the value of 77 in the random bond-disordered three- 
state Potts ferromagnet does not change with the strength of disorder, while both v and 7 
vary continuously. Our results disprove so-called semi-rigorous inequality, but satisfies the 
concept of weak universality. The features found in the model seem to be generic in any 
disordered ferromagnets irrespective of the value of a p , contrary to the conventional wisdom. 
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Figure Captions 

Figure(l): ln(x^) as a function of ln(L). The slope of a straight line would correspond to 
2 — rj. The dotted line represents the case of the exact value of rj of the undisordered system. 
Figure(2): Xl/Xoo as a function of at certain fixed temperatures. For each J', xl 

becomes equivalent to its bulk value within the statistical errors under the condition L/^ L > 
5. Note that this condition is independent of temperature according to the standard theory 
of finite size scaling, that is, e.g., Xl(P)/Xoo(P) = q x (^L(P)/L). 

Figure(3): ln(x(t)) as a function of | ln(t)|. The slope of a straight line is equivalent to the 
value of 7. Each dotted line represents the best approximation of the data obtained from 
the x 2 fit- R is clear that the value of the slope increases with decreasing J' for J' < 0.9, 
and that for each J' it does not vary with t. 
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